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ABSTRACT 

We present ISO LWS observations of the proto-planetary nebula CRL 618, 
a star evolving very fast to the planetary nebula stage. In addition to the lines 
of i^CO , i^CO , HON and HNC, we report on the detection of H2O and OH 
emission together with the fine structure hues of [01] at 63 and 145 fj.m. The 
abundance of the latter three species relative to ^^CO are 4 10~^, 8 10""^ and 4.5 
(approximate value) in the regions where they are produced. We suggest that 
0-bearing species other than CO are produced in the innermost region of the 
circumstellar envelope. The UV photons from the central star photodissociate 
most of the molecular species produced in the AGB phase and allow a chemistry 
dominated by standard ion-neutral reactions. Not only allow these reactions 
the formation of 0-bearing species, but they also modify the abundances of 
C-rich molecules like HCN and HNC for which we found an abundance ratio of 
~ 1, much lower than in AGB stars. The molecular abundances in the different 
regions of the circumstellar envelope have been derived from radiative transfer 
models and our knowledge of its physical structure. 

Subject headings: infrared: stars — ISM: molecules — line: 
identifications — planetary nebulae: individual (CRL618) — stars: 
abundances — stars: carbon 
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1. Introduction 

CRL 618 is one of the few clear examples of an AGB star in the transition phase to 
the Planetary Nebula stage. It has a compact HII region created by a hot central star 
(Westbrook 1976, Kwok & Feldman 1981). This object is observed as a bipolar nebula 
at optical, radio and infrared wavelengths (Carsenty & Solf 1982, Bujarrabal et al. 1988, 
and Hora et al. 1996). The expansion velocity of the envelope is around 20 kms"'^, but 
CO observations show the presence of a high-velocity outflow with velocities up to 300 
kms~^(Cernicharo et al. 1989). High-velocity emission in the B.2 v = 1 — S(l) line is 
also detected (Burton & Geballe 1986). The high velocity wind and the UV photons from 
the star perturb the circumstellar envelope (CSE) producing shocks and photodissociation 
regions (PDRs) which modify the physical and chemical conditions of the gas (Cernicharo 
et al. 1989, and Neri et al. 1992). 

We present in this Letter an ISO (Kessler et al. 1996) Long Wavelength Spectrometer 
(LWS, Clegg et al. 1996, Swinyard et al. 1996) spectrum between 43 and 197 /im. In 
addition to the high-J lines of CO, HCN and HNC, we have discovered, for the flrst time 
in a C-rich star, several lines of water vapor, OH and the flne structure lines of [01]. We 
discuss the origin of these 0-bearing species and we present a model that accounts for the 
observed emission. 

2. Observations and Results 

The LWS data were taken on revolution 688. The data have been processed following 
pipeline number 7. We have used ISAP^ to remove glitches and fringes. The rotational lines 
(see Fig. 1 and 2) of ^^CO (J=14-13 to J=41-40), of ^^CO (J=14-13 to J=19-18), HCN, 
HNC, [01] (at 63.170 and 145.526 /xm) and C+ (^^3/2 pi/2 transition at 157.74 /im) are 
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observed. Several lines of OH and H2O are also detected (see Fig. 1 and 2). 

The far-infrared spectrum is essentially dominated by the CO lines. All the other 
species have much lower intensities contrary to what is found in the AGB star IRC-l-10216 
where the HCN hues (from the ground and vibrationally excited states) are as strong as 
CO. The main difference is due to the physical structure of the CSE of CRL 618 which has 
a central hole in molecular species filled by a bright HII region. The LWS data also show 
several weak hues that remain unidentified. We have searched for CH+, NH, CH, and other 
light species. But unlike the case of NGC7027 (Cernicharo et al. 1997) where the PDR has 
a large spatial size, the PDR of CRL 618 is still close to the central object and has a larger 
dilution in the LWS beam, which makes the fluxes of these expected species very weak. 

The full spectrum of CRL 618 between 2 and 200 /im has been discussed by Cernicharo 
et al. (1999a et 1999b). In these data, in addition to the species detected in this Letter, they 
have also reported the presence of several long polyacetylenic chains and the bending modes 
of HC3N and HC5N. However, in the LWS part of the spectrum no pure rotational lines of 
HC3N have been found. Nevertheless, the i/7 band of HC3N at 44 /xm is observed and these 
data will be published elsewhere. All the observed absorption bands in the mid-infrared 
arise from gas at ~ 250 K and densities larger than 10^ cm~^ (Cernicharo et al. 1999b). 

3. Discussion 

Previous works on CRL 618 underlined the presence of a central torus associated 
with a Photo-Dissociated Region, an extended AGB remnant envelope and high velocity 

^The ISO Spectral Analysis Package (ISAP) is a joint development by the LWS and SWS 
Instrument Teams and Data Centers. Contributing institutes are CESR, IAS, IPAC, MPE, 
RAL and SRON. 
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wind regions (HVW) as bipolar lobes (Weintraub et al. 1998, Burton & Geballe 1986, 
Neri et al. 1992, Carsenty & Solf 1982), whose physical parameters, like the velocity of 
the gas or the angular size, were investigated (Martin-Pintado et al. 1995, Meixner et 
al. 1998, Hajian et al. 1996). In order to reproduce the emission, we have to introduce 
several spatial components, in terms of kinematics, temperature and density based on 
the knowledge we have of this object. In Fig. 3 we show the adopted geometry: (i) a 
central torus with the PDR; (ii) the extended AGB remnant envelope; (iii) the lobes 
(bipolar jets, HVW). We tried to reproduce the LWS observations by simulating a spectrum 
(12C0 +^^C0 +HCN+H2O +OH+[OI]). We first tried to model the ^^CO intensities. 
Then we applied the same parameters to the other molecules allowing to vary only these 
column densities (see parameters in Table 1). We used a simple LVG-model and obtained 
satisfactory fits for several molecules (see Fig. 1). 

The torus is assumed to be composed of a large neutral shell of 1.1"-1.5" size, 
temperatures from 250 to 800 K, density of 5 10^ cm~^, and an expansion velocity of 20 
kms~^. All the detected molecules are present in the 250 K area where most of the CO long 
wavelength emission comes from. In this region the ^^CO column density (TV) is 10^^ cm~^ 
and the ^^CO , H2O , HON, HNC and OH abundances relative to ^^CO are 1/20, 1/25, 
1/1000, 1/1000 and 1/1250 respectively. We adopt an ortho-to-para ratio for H2O around 
3. A warmer region in the torus is needed to reproduce the ^^CO and ^^CO emissions 
at shorter wavelengths (size ~ 0.6",N= 1 10^^ cm~^, T=1000 K) and the atomic oxygen 
emission. As we reach this inner edge of the PDR, only the atomic oxygen seems to be 
present, whereas all molecular species are photodissociated: N ([01]) =4.5 10^^ cm~^, thus 
an abundance of 4.5 relative to CO. But this value must be considered as very approximate 
due to the uncertainties in the size of this region. For the inner radius of the torus, which 
should be in contact with the HII region, we calculate at a distance of 0.3" (10^^ cm) a UV 
flux of 20 000 Go (in units of the local IS radiation flux); this UV illumination is around 
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200 in the lobes (Latter et al.l992). 

The contribution from the bipolar lobes corresponds to the high velocity gas (200 
kms~^, Cernicharo et al. 1989). Previous observations of CO line wings (Cernicharo et 
al. 1989) and H2 lines (Hcrpin et al. in preparation) indicate high velocity and dense 
shock-heated molecular gas in CRL 618. The contribution of such a region of dense warm 
shocked gas could dominate the CO high rotational transitions (e.g., HoUenbach & McKee 
1989). The high temperature component at high velocity is necessary to produce part of 
the emission at short wavelengths, since a high temperature component arising from the 
torus produces too weak intensities. We assume that the bipolar lobes in CO have a size 
between 1.5" and 1.7" with Arranging from 2 10^^ to 5 10^^ cm~^, temperatures of 200 K 
and 1000 K, and density around 10^ cm~^. A contribution from this HVW region for H2O , 
OH and [01] line emission is likely because the molecules in the AGB remnant must be 
dissociated in the shocks, and in the molecular recombination some 0-bearing species could 
be found like in the inner part of the torus. However, due to the lack of velocity resolution, 
we can not discriminate both components. 

The contribution of the extended AGB remnant envelope is as a colder zone (100 K) 
for ^^CO and ^^CO . The density is around 5 10^ cm^^ and can be compared to what was 
found by Yamamura et al. (1994) with a mass loss rate of 10~^ Mq/jt for CRL 618, and 
Neri et al. (1992) (10^ - 10^ cm~^ on the outer edge of the lobes). 

For the HVW region and the AGB remnant envelope, we adopt a temperature profile 
derived from the AGB phase case (Langer & Watson 1984, Herpin 1998), going from 1000 
K on the origin of the conical structure (corresponding to the temperature at 6 — 7 10^^ cm 
for a Mira at 3000 K) to 100 K on the exterior, equal to the AGB temperature for dust 
in the outer part of the envelope (distance > 10^^ cm) and equal to the value derived by 
Yamamura et al. (1994). 
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The model fits well all the CO lines, with small discrepancies at short wavelengths 
mainly due to the presence of different molecular contributions, to a lower signal-to-noise 
ratio, and to the fact that some CO lines are mixed with H2O and OH lines. Nevertheless, 
all the observed line intensities from J=14-13 up to J=41-40 are reproduced. Our model also 
fits the observed emissions J=l-0, 2-1 (Cernicharo et al. 1989), 3-2 (Gammie et al. 1989) 
and 6-5 (unpublished data). A good fit is also obtained for the atomic oxygen lines at 63 
/im and 145 //m. 

Although the ^^CO emission is weak, the full coverage of the LWS spectrum allows us 
to derive a [^^CO /^^CO ] ratio of 20, lower than the common interstellar value of 40-45, 
but equal to the detected ratio (20) of Kahane et al. (1992) in this object. 

Emission from the HCN wings (Neri et al. 1992) is produced by the strong impact of a 
primary outflow on a few localized dense clumps in the inner envelope. According to these 
authors, 70% of the HCN line emission comes from a very compact, spherically symmetric 
region, and is mainly coUisionally excited and strongly affected by self-absorption. They 
found for the flow (our HVW region) n{H2) = few 10^ - 10^ cm'^ N{CO)^ 2 10^^ cm'^ 
and N{HCN) = few 10^^ cm-2, which leads to [HCN/CO] around 10"^ Our high-velocity 
HCN component (see Table 1) comes from the lobes (HVW) (A^ = 2 10^^ cm'^, T= 1000 
K). We find here a ratio [HCN/CO] of 10"^ equal to the ratio obtained by Neri et al. (1992). 
When one considers an abundance of CO of a few 10~^, this leads to a HCN abundance of a 
few 10~^ for our data ([HCN/HNC] ~ 10). This is the expected abundance ratio in a C-rich 
star in the AGB phase (Cernicharo et al. 1996), and proves that this HCN material has 
not been processed by the UV photons from the central hot star. Furthermore, Deguchi et 
al. (1986) found in the disk surrounding the source an abundance of HCN of order of 2 10~^ 
per H2 molecule, much less than Neri's value. For this region (sec Table 1) we have found 
a ratio [HCN/CO] of 10~^ which leads to an abundance of a few 10~^. We can conclude 
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that in this region HCN molecules are already efficiently dissociated. The HCN and HNC 
abundances are similar, which implies a high degree of reprocessing of the molecular gas 
(x(HCN) » x(HNC) in standard C-rich AGB stars). 

Below 110 fxm our models predict that the H2O and OH lines will be in emission. 
However, the data seem to indicate that some of these lines could be in absorption or 
are too weak. Taking into account the reduced sensitivity at these wavelengths and the 
complex structure of the region, a more detailed model including radiative transfer effects 
for H2O and OH, similar to those discussed by Gonzalez-Alfonso & Cernicharo (1999), will 
need a better knowledge of the physical and geometrical structure of the innermost region 
of the envelope. 

The detection of the fine structure lines of [01] at 63 and 145 /im gives important 
insight into the chemical evolution of CRL 618. This indicates that the reservoir of CO 
molecules formed in the AGB phase of CRL 618 has started to be reprocessed through UV 
photons and shocks. In CRL 618 Cernicharo et al. (1989) also found H2CO and Bujarrabal 
et al. (1988) found HCO+. These molecules are not formed efficiently in the AGB phase 
of a C-rich star. How are they formed in CRL 618 ? As soon as some atomic oxygen 
is released into the gas phase from the photodissociation of CO, standard ion- neutral 
reactions can form quickly H3O"'" and HCO"'". From the former, and through electronic 
recombination, H2O and OH are efficiently produced. We expect these molecules to arise 
from the innermost region of the envelope where the UV flux is largest. In our models, 
the column densities in that region are 10^^, 5 10^^ and 8 10^^ cm~^ for ^^CO , H2O and 
OH respectively. That means that the abundances of H2O and OH relative to CO are 
4 10"^ and 8 10"^ respectively. These abundances are much lower that those found in the 
surrounding media of galactic HII regions (Cernicharo et al. 1999a) and suggest that only 
a small fraction of the UV photons from the central star are reaching the region where 
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H2O and OH are formed. In fact, we could expect a dusty region in the torus around 
the central HII region that protects partially the molecules formed in the AGB phase. 
Obviously, in the PDR itself all molecules could be destroyed and most of the [01] dominant 
coolant in high density PDRs emission is probably coming from that region. 

The model derived from the LWS spectra leads to two global regions of molecular 
emission. One in the torus with several molecular species (^^CO , ^^CO , HON, HNC, 
H2O and OH) at moderate temperature. Another in the lobes dominated by CO and HCN 
molecules. The detected atomic oxygen line is strong and comes from the innermost warm 
region, the inner edge of the PDR, what proves that CRL 618 is already quite well on its 
way to the PN stage and that these objects, even in the C-rich case, can produce efficiently 
0-bearing molecular species. 

Acknowledgements We thank Spanish DOES and CICYT for funding support for this 
research under grants PB96-0883 and ESP98-1351E. 
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Table 1: Table of the column densities (in cm~^) for all molecules present in each region of 
the object as depicted in Figure 3 (see text for the adopted parameters). 

Molecules ^^CO ^^CO HCN HNC 0-H2O OH 0[I] 

TORUS (including PDR) 

T=250 K,0 = 1.5" 

V = 20 kms-i, n(if2) =5 10^ cm-^ 
T=800 K, 6* = 1.1" 

V = 20 kms-i, n{H2) =5 lO'^ cm-^ 
T=1000 K, 6' = 0.6" 

V = 20 kms-i, n(if2) =7 lO'^ cm-^ 
HVW 

T=200 K, 6' = 1.7" 

V = 50 kms-i, n(i?2) = lO'^ cm-^ 
T=1000 K,e = 1.5" 

V = 200 kms^i, n{H2) = 10^ cm-^ 
AGB remnant 

T=100 K, 6' = 5" 7 10^^ 3.5 10^^ 7 10^^ 7 lO^^ 

V = 20 kms-i, n(i?2) = 5 10^ cm-^ 



lO^^ 5 IQi^ 10^6 10^6 3 IQi^ 8 10^5 
6 10^^ 3 IQis 

1 10^9 5 IQi'^ 4.5 10^^ 



5 IOI8 2.5 10^^ 3 IQi^ 2 10^7 
2 10^6 1 IQis 2 IQis 1014 



Figure Captions 



Fig. 1. — Continuum subtracted spectra of CRL 618 (LW detectors). The result of our 
model is shown by the continuous hne. The lines of CO, ^^CO , HCN, H2O and OH are 
indicated by arrows while those of HNC are indicated by vertical lines (from J=22-21 at 
150.627 fim. to J=17-16 at 194.759 /xm). The C+ transition is not included in our models; 
the plot indicates a gaussian fit to this feature. 

Fig. 2. — Important detected lines in CRL 618 (continuum subtracted spectra). 

Fig. 3. — Simple view of the molecular and emission structure as seen by the model. HVW 
is for the High Velocity Wind region. 
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